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ABSTRACT 

Recently, high-resolution Chandra observations revealed the existence of very sharp features in the X- 
ray surface brightness and temperature maps of several clusters (Vikhlinin et al., 2001). These features, 
called cold fronts are characterized by an increase in surface brightness by a factor > 2 over 10-50 kpc, 
accompanied by a drop in temperature of a similar magnitude. The existence of such sharp gradients 
can be used to put interesting constraints on the physics of the intracluster medium (ICM), if their 
mechanism and longevity are well understood. Here, we present results of a search for cold fronts in 
high-resolution simulations of galaxy clusters in cold dark matter (CDM) models. We show that sharp 
gradients with properties similar to those of observed cold fronts naturally arise in cluster mergers when 
the shocks heat gas surrounding the merging sub-cluster, while its dense core remains relatively cold. The 
compression induced by supersonic motions and shock heating during the merger enhance the amplitude 
of gas density and temperature gradients across the front. Our results indicate that cold fronts are non- 
equilibrium transient phenomena and can be observed for a period of less than a billion years. We show 
that the velocity and density fields of gas surrounding the cold front can be very irregular which would 
complicate analyses aiming to put constraints on the physical conditions of the intracluster medium in 
the vicinity of the front. 

Subject headings: cosmology: theory - intergalactic medium - methods: numerical - galaxies: clusters: 
general - instabilities-turbulence-X-rays: galaxies: clusters 



1. INTRODUCTION 

Recently, high-resolution Chandra observations revealed 
the existence of very sharp features in the X-ray sur- 
face brightness and temperature maps of several clus- 
ters, including A2142 (Markevitch et al., 2000) and A3667 
(Vikhlinin et al., 2001). The sharp gradients in the X-ray 
surface brightness were first found in the ROSAT images 
of A2142 and A3667 along with signatures of substruc- 
tures (i.e., several small groups and merging components 
embedded in a larger cluster). Due to lack of spatially 
resolved temperature maps, these features were initially 
interpreted as shock fronts arising during cluster mergers 
(Markevitch et al., 1999). This interpretation turned out 
to be incorrect when detailed X-ray surface brightness and 
temperature maps were obtained with the Chandra satel- 
lite. The maps showed that the increase in gas density 
and X-ray surface brightness by a factor of > 2 across 
the sharp features is accompanied by a decrease in tem- 
perature of gas of a similar magnitude, the behavior op- 
posite to that expected across shock fronts. These sharp 
features were therefore named ''cold fronts". Markevitch 
et al. (2000) and Vikhlinin et al. (2001) interpreted these 
fronts as a boundary of "the dense subcluster core that 
has survived a merger and ram pressure stripping by the 
surrounding shock- heated gas." 

The existence of the sharp temperature gradients can be 
used to put interesting constraints on the conditions in the 
intracluster medium (ICM) in the vicinity of the front. In 
particular, the width of the observed front (< 5 — 10 kpc) 
is several times smaller than the Coulomb mean free path 
for electrons in the ICM. This indicates that the thermal 
conduction must be suppressed, at least across the front. 
In the case of A2142, Ettori & Fabian (2000) find that the 



thermal conductivity has to be reduced by several orders of 
magnitude from the classical Spitzer value near the front. 
The observed extent, ~ 500 kpc, of the sharp boundary 
may also indicate that the Kelvin-Helmholtz instabilities 
are partially suppressed, most likely by magnetic fields 
parallel to the boundary. The instabilities, expected to 
arise on the boundary if the hot gas flows along the front, 
would disrupt and widen the boundary. Using a simple dy- 
namical model that approximates the gas cloud as a dense 
spherical cloud, Vikhlinin et al. (2001) find that the mag- 
netic field of 7-16/iG parallel to the front must be present 
at the boundary to maintain its stability. The detailed 
studies of cold fronts can, thus, provide new detailed in- 
sights into the physics of the ICM, including the efficiency 
of energy transport and the magnetic field strength in the 
ICM. The reliability of such models, however, depends on 
understanding the dynamics of gas motions during cluster 
mergers and, especially, in the vicinity of the cold fronts. 
The constraints on the transport mechanisms require also 
that we understand how long the cold fronts can survive 
dynamically. 

In this paper we present results of the search and anal- 
ysis of features similar to the observed cold fronts in very 
high-resolution cosmological N-body+gasdynamics simu- 
lations of clusters forming in Cold Dark Matter (CDM) 
models. The main goals of this paper are 1) to test the 
interpretation of the cold fronts as boundaries between hot 
shock-heated gas and dense cold gas of the merging sub- 
clump; 2) identify the situations and mechanisms that can 
produce cold fronts and 3) analyze the dynamics and den- 
sity fields of gas and dark matter in the vicinity of the 
front. 

Although numerous extensive theoretical analyses of 
cluster mergers have been done in the past, the typical 
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spatial resolution in studies of cluster formation in the 
full cosmological context is too low to match the superb 
spatial resolution of the Chandra images. Eulerian gasdy- 
namics codes based on the modern high-resolution shock 
capturing techniques are capable of resolving sharp density 
and temperature gradients of arbitrary amplitude within 
1 — 2 grid cells and are thus well suited for studies of 
shocks and cold fronts in cluster mergers. However, high- 
dynamic range required in self- consistent simulations of 
cluster formation limited use of the Eulerian codes to con- 
trolled merger experiments (e.g., Roettiger et al., 1997; 
Quilis et al., 1998; Ricker & Sarazin, 2001). The results of 
these simulations were widely used for physical interpre- 
tation of the new high-resolution X-ray observations of 
merging clusters. More recently, the advent of cosmologi- 
cal codes using the adaptive mesh refinement (AMR) (e.g., 
Bryan & Norman, 1997; Kravtsov et al., 2002; Teyssier, 
2002) allowed to achieve the dynamic range of ^ 10^ — 10^ 
with Eulerian gasdynamics algorithms, thus making pos- 
sible high-resolution self-consistent cluster simulations in 
a realistic cosmological setting. As we will show below, 
these simulations can resolve sharp gradients in density 
and temperature fields of the ICM on scales of ~ 10 kpc, 
approaching the typical resolution of Chandra images. 

The paper is organized as follows. In §2, we discuss 
the cluster simulations used in our analysis. In §3, we 
discuss and illustrate the detailed structural evolution of 
the cluster gas and surrounding filaments during cluster 
mergers and accretion of matter. In § 4 we present cold 
fronts identified in numerical simulations and discuss their 
origin, structure, dynamical properties, and detectability. 
We discuss our results and summarize our conclusions in 
§5. 

2. NUMERICAL SIMULATIONS 

In this study, we analyze the high-resolution cluster sim- 
ulations performed using the Adaptive Refinement Tree 
(ART) A^-body+gasdynamics code (Kravtsov, 1999; Kravtsov 
et al., 2002). ART is an Eulerian code designed to achieve 
high spatial resolution by adaptively refining regions of 
interest, such as high-density regions or regions of steep 
gradients in gas properties, and has good shock-capturing 
characteristics. The code can capture discontinuities in 
gas properties within ~ 1 — 2 grid cells, which makes it 
well-suited for studying sharp features such as shock and 
cold fronts. At the same time, the adaptive refinement 
in space and (non-adaptive) refinement in mass (Klypin 
et al., 2001) allows us to reach the high dynamic range 
required for high-resolution self- consistent simulations of 
cluster evolution in a cosmological setting. 

We note that the effects of magnetic fields, gas cool- 
ing, stellar feedback, and thermal conduction are not in- 
cluded in these simulations. This is the first theoretical 
study of cold fronts and it is therefore interesting to see 
whether these features can be produced in the case of 
purely adiabatic gasdynamics. As we show below, cold 
fronts arise at the interfaces between relative low density 
hot medium and dense cold gas of merging clumps. These 
regions should not be significantly affected by gas cooling, 
which is expected to have the largest effect in cluster cores. 
The absence of thermal conduction is acceptable because 
our primary goal is to study the dynamical origin of cold 
fronts, rather than their precise properties. The very ex- 



istence of cold fronts indicates that thermal conduction is 
suppressed in their vicinity. Other theoretical arguments 
(e.g., Loeb, 2002) constrain the thermal conduction coeffi- 
cient to be significantly smaller than the classical Spitzer 
value over most of the ICM, except perhaps in the cluster 
cores (Voigt et al., 2002). Our simulations, in fact, implic- 
itly assume that the Coulomb mean free path is smaller 
than the size of computational cells (this is one of the 
main assumption of gasdynamics), and hence there is no 
thermal conduction. 

We analyze simulations of two clusters of intermediate 
mass in the ACDM and standard CDM models. The first 
(ACDM) cluster was simulated in the "concordance" flat 
ACDM model: Qrn = I - = 0.3, = 0.043, h = 0.7 
and <J8 = 0.9, where the Hubble constant is defined as 
lOO/i km s~^ Mpc~^, and as is the power spectrum nor- 
malization on Sh~^ Mpc scale. The simulation used 128^ 
uniform grid and 7 levels of mesh refinement in the com- 
putational box of SOh~^ Mpc, which corresponds to the 
dynamic range of 128 x 2^ = 16, 384 and peak resolution 
of 80/16,384 ^ bh-^ kpc. Only the region of - lOh'^ 
Mpc around the cluster was adaptively refined, the rest 
of the volume was followed on the uniform 128^ grid. Sim- 
ilarly, the simulation used initial conditions sampled with 
different particle masses (see Klypin et al., 2001). The 
mass resolution (dark matter particle mass) in the region 
around the cluster was 2.7 x 10^h~^ M©, while other re- 
gions were simulated with lower mass resolution (the total 
of three particle species, different in mass by a factor of 
eight, were used). The region of the cluster was selected 
for re-simulation at higher resolution from a lower resolu- 
tion AT-body run. 

The second (SCDM) cluster is simulated from the initial 
conditions used in the "Santa Barbara cluster comparison 
project" (Frenk et al., 1999). The cosmological parame- 
ters assumed were those of the standard CDM cosmology: 
Qru 1.0, 1^6 0.1, h 0.5 and as 0.9. The simu- 
lation used 128^ uniform grid and 6 levels of refinement, 
which corresponds to the peak resolution of ~ 4/i~^ kpc 
in the 32h~^ Mpc computational box. Here, again we 
used a resampling technique to set up the initial conditions 
for multiple particle species. The lowest particle mass is 
4.9 X 10^ Mq and corresponds to 256^ effective num- 
ber of particles. The Lagrangian region corresponding to 
the spherical volume of the radius R = 5i?vir, where i^vir 
is the virial radius of the cluster at z = 0, was simulated 
with the highest mass resolution and was adaptively re- 
fined. The SCDM cluster simulation was previously used 
to compare performance of the ART code to other existing 
cosmological gasdynamics codes in a realistic cosmological 
setting (Kravtsov et al., 2002), where good agreement was 
found in both cluster properties and its radial profiles. 

The SCDM and ACDM clusters had virial masses of 
M200 = 6.0 X lO^'^h-^ Mo and M340 = 2.4 x lO^'^h'^ M© 
and virial radii of R200 = 1.35/i~^ Mpc and R340 = 
1.26h~^ Mpc, respectively. Here subscripts indicate the 
cosmology-dependent virial overdensity with respect to the 
mean density of the universe. Both clusters were therefore 
simulated with about one million dark matter particles 
within the virial radius. In addition to 128^ cells of the 
uniform grid, approximately 10^ mesh cells were generated 
in the process of adaptive refinement in the regions sur- 
rounding clusters. The simulations were thus designed to 
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Fig. 1. — The maps of projected density (top) and entropy (defined as Tn ^Z"^, bottom) of a simulated ACDM cluster at four different 
redshifts in a 60h~^ kpc slice centered on the central density peak. The maps are color-coded on a log^^Q scale in units of cm^ (density) and 
keV cm^ (entropy). The size of the region shown is 8h~^Mpc (10% of the entire simulation volume). The entropy maps reveal a very complex 
entropy distribution of the gas. Both the filaments and the forming cluster are surrounded by strong accretion shocks. Note, however, that 
the accretion shock around the cluster is very aspherical and does not penetrate into the filament; relatively low-entropy gas accreting onto 
cluster along the filament does not pass through the strong virial shocks and can be traced all the way to the central 0.5h~^ Mpc of the 
cluster. 



concentrate the bulk of the computational effort on the 
individual clusters, which allowed the simulations to reach 
the high dynamic range critical for this study. The spatial 
resolution of these simulations allows us to resolve sharp 
gradients over scales < 20h~^ kpc within the virial radius 
of clusters. This is comparable to, albeit still somewhat 
worse than, the superb angular resolution of the Chandra 
satellite. 

3. CLUSTER EVOLUTION 

Let us start the discussion of cold fronts with a gen- 
eral description of cluster evolution. Figure 1 shows maps 
of projected gas density (top) and entropy (defines as 
Tn~^/^, bottom) of the simulated ACDM cluster at four 
different redshifts in a 60h~^ kpc slice centered on the 
central gas density peak. The maps are color-coded on 
a log;LO scale in units of cm^ (density) and keV cm^ (en- 
tropy). The size of the region shown is 8/i~^Mpc (10% of 
the entire simulation volume). The figure illustrates com- 
plex dynamical processes accompanying cluster formation: 
accretion of clumps and diffuse gas along filaments, strong 
accretion and virial shocks both around cluster and sur- 
rounding filaments, weaker merger shocks within the virial 
shock of the cluster, and the complicated fiow pattern of 
the ICM gas revealed by the entropy map. 

The evolution of intergalactic gas is driven by accretion 
and mergers between clumps. Between z = 1 and z = 0.5, 
the main cluster undergoes a nearly equal mass merger. 
Such mergers are spectacular events involving the kinetic 
energy as large as ~ 10^^ ergs, the most energetic events 
since the Big Bang. More common are smaller mergers 
and accretion of groups along large-scale filaments. At 



z ~ 1 — 3, highly supersonic accretion (Mach numbers of 
> 100) of pristine gas from low-density regions leads to 
formation of strong shocks located roughly at the virial 
radii of the two merging proto-clusters. Note that shocks 
of similar strength also surround the filament along which 
these clusters move. Note also that virial shocks do not 
penetrate into the filament itself and it is in fact difficult 
to separate the virial shocks of clusters from the accretion 
shocks around the filament. 

During the evolution from z = 1 to z = 0, as the mass 
of the main cluster grows, the radius and entropy gradient 
of the virial shocks steadily increase. Similar increases can 
be observed for the accretion shocks around the filament. 
However, even at the present epoch the gas fiowing along 
the filament reaches the central regions of cluster without 
passing through a strong virial shock. This is clearly illus- 
trated in the entropy maps in which relatively low-entropy 
gas flowing along the filament can be traced all the way to 
the central 0.5h~^ Mpc of the cluster. The virial shocks 
instead propagate along directions of the steepest pressure 
gradient into the low-density voids. The flow of fllamen- 
tary gas, with entropy already increased by the accretion 
shocks (temperature of several million degrees K), has con- 
siderably smaller Mach numbers (typically ~ 1 — 10). Its 
motion is therefore slightly supersonic at best, compared 
to the highly supersonic motion of the yet unshocked gas 
flowing from voids. When this gas reaches the cluster core 
it accelerates and generates random, weak shocks in which 
it dissipates its kinetic energy. The flows from different 
fllaments and directions result in a complicated flow pat- 
terns. We will discuss these motions in relation to the cold 
fronts below (see also Figure 5). 
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Fig. 2. — The X-ray surface brightness (top) and emission- weighted temperature (bottom) of the simulated ACDM cluster at four different 
redshifts. The maps are color-coded on a log^o scale in units of erg s"-*^ cm~^ arcmin"^ (surface brightness) and keV (temperature). Both 
surface brightness and temperature are calculated in the 0.5-2 keV band. The size of the volume shown is 2h~^ Mpc. Note the cold front 
associated with the merging sub-cluster behind the merger shock front to the left of two sub-clusters in the z = 0.43 panel. Note also the 
cold front associated with the merging subclump at z = to the left of the cluster center and the adjacent compressed region of enhanced 
temperature. The subclump is trailed by a relatively cold (~l-2 keV) intergalactic gas accreted along a filament. 



4. COLD FRONTS 

The main question we want to answer in this study is 
whether cosmological simulations of galaxy clusters form- 
ing in hierarchical models can produce features similar to 
the "cold fronts" observed by the Chandra satellite. If 
so, simulations can then be used to shed light on the ori- 
gin and properties of these features. In this section, we 
analyze the high-resolution cluster simulations described 
above to search for sharp gradients in temperature and 
density structure of clusters. 

4.1. X-ray surface brightness and temperature maps 

To perform the search and analysis of cold fronts, we 
generate the X-ray surface brightness and emission- weighted 
temperature maps of the simulated clusters. The maps 
were generated from uniform density and temperature 
grids centered on the minimum of cluster potential. The 
density and temperature grids, in turn, were constructed 
using interpolation from the original mesh refinement 
structure. The X-ray emissivity of gas cells was calcu- 
lated in the 0.5-2 keV energy band using the Raymond & 
Smith (1977) plasma emission model and assuming a uni- 
form metallicity of 0.3 solar. Figure 2 shows the X-ray sur- 
face brightness (top) and emission-weighted temperature 
(bottom) maps of the simulated ACDM cluster at four dif- 
ferent redshifts. The maps are color-coded on a log^g scale 
in units of erg s~^ cm~^ arcmin"^ (surface brightness) and 
keV (temperature). The region shown in the figure is 2h~^ 
Mpc; the pixel size is 2h-^ Mpc/256 ^ Sh'^ kpc. The 
size of the pixels were chosen to be similar to the effective 
resolution of the simulation within the cluster. 

These maps show that complex density and tempera- 
ture structure of the intracluster gas, illustrated by the 



thin density and entropy slices in Figure 1, is present in 
the observables averaged over the line-of-sight. At z = 1, 
two clusters are approaching each other along a massive 
filament. These clusters undergo a nearly equal-mass, 
slightly off- axis merger at z = 0.54. Compression of gas 
between two merging clumps significantly enhances the 
temperature in the large band visible in z = 0.54 panel. 
By z = 0.43, the clusters have passed each other once 
and their supersonic motion results in two large-scale bow 
shocks propagating in opposite directions along the merger 
axis. These rather strong shocks are clearly seen in the 
temperature map of z = 0.43 epoch as regions of enhanced 
temperature with very sharp boundaries. At z = 0, the 
cluster formed by the z ^ 0.5 merger is relatively relaxed 
but continues to accrete matter and undergo minor merg- 
ers (note a small merging sub-clump to the left of the 
cluster center). 

4.2. Origin, properties and detectability 

The maps in Figure 2 show that a variety of very sharp 
features in temperature and surface brightness arises dur- 
ing minor and, especially, major mergers. Some of these 
features, indeed, exhibit properties qualitatively similar to 
the observed cold fronts. Below, we present two cases of 
cold fronts found in the CDM simulations. 

4.2.1. Large-scale cold front during a major merger 

A spectacular cold front arises in a major merger of 
the ACDM cluster dit z ^ 0.4 - 0.5. This cold front is 
a low temperature region embedded in the high tempera- 
ture gas just behind the left bow shock in the z = 0.43 
panel of Fig. 2. The colder gas belongs to one of the 
merging clusters. The dark matter halo of this cluster 
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Fig. 4. — The mock Chandra photon count image of the cold front 
visible in the z = 0.43 panel of Figure 2. The image was constructed 
assuming 50 ksec exposure in the 0.5—4 keV band, cluster redshift of 
z = 0.05, and a flat background level of 4.3 x 10 ~^ cuts s"-*^ pixel" 
The flgure shows the full 16' x 16' ACIS-I fleld of view binned to 2'' 
pixels. The color represents photon counts per pixel in logio scale. 



Fig. 3. — The X-ray surface brightness (top), emission- weighted 
temperature (middle) and emission-weighted pressure (bottom) pro- 
flies across cold fronts identifled in z=0.43 (left) and z=0 (right) 
panels in the Figure 2. The distance in kpc is measured relative 
to the cold fronts shown by dashed lines. In both cases, a jump in 
surface brightness is accompanied by a drop in temperature over the 
scale of ~ 30 kpc, while pressure changes relatively smoothly across 
these features. The simulated proflles and images reproduce all the 
main features of the observed cold fronts. The proflles in the left 
panels also show the presence of a bow shock ^ 200 kpc from the 
cold front. 

is in the process of disruption and the cluster gas ap- 
pears to be "sloshed out" of its potential well (see Fig. 5). 
Rapid motion of the gas inside the freshly shocked gas 
then results in a sharp boundary between cold unshocked 
gas of the cluster and the hot shocked surrounding gas. 
The left panels of Figure 3 show the X-ray surface bright- 
ness (top), emission- weighted temperature (middle) and 
emission-weighted pressure (bottom) profiles across this 
temperature gradient in the z = 0.43 maps. The pro- 
files were constructed by drawing a straight line across the 
front and selecting values of map pixels nearest to the line 
in equally spaced intervals. The distance in kpc is mea- 
sured relative to the sharp temperature gradient indicated 
by the dashed line. The profiles show that a drop of tem- 
perature over the scale of 20 — 30h~^ kpc is accompanied 
by a jump in surface brightness by a factor of ^ 2, the 
behavior is very similar to the profiles of observed cold 
fronts. Despite jumps in gas density and temperature, 
pressure changes relatively smoothly across this bound- 
ary. The smooth change of pressure was also deduced for 
observed cold fronts. Note that the profiles show the pres- 
ence of a bow shock ^ 200 kpc from the cold front. A 
similar bow shock was also observed in A3667 ahead of 
the cold front, although its detection was marginal. 

The overall structure of the front also exhibits qual- 
itatively similar behavior to those observed in A3667 



(Vikhlinin et al., 2001). To illustrate its appearance in 
a realistic observations, we constructed the mock photon 
count Chandra X-ray map of the region around the sim- 
ulated cold front. The map is shown in Figure 4. We 
assumed a 50 ksec exposure in the 0.5 — 4 keV band, clus- 
ter redshift of z = 0.05, and a flat background level of 
4.3 X 10-^ cuts s-i pixel-^ The fi gure shows the full 
16' X 16' ACIS-I field of view binned to 2'' pixels. The color 
represents photon counts per pixel in logio scale. The fig- 
ure shows that the cold front is detectable by Chandra and 
appears remarkably sharp in the photon count map. The 
boundary of the front spans ~ 0.5h~^ Mpc, comparable 
to the size of the observed cold front in A3667. At larger 
angles, the sharp discontinuity widens and then gradually 
disappears, the onset of Kelvin-Helmholtz instability also 
seen in A3667; however, the appearance of the simulated 
cold front seems to be more irregular than that of the ob- 
served front (see also Fig. 5 and §5 for more discussions). 

To examine the processes that produce and then destroy 
this cold front, we created a movie of the temperature, 
density, and X-ray surface brightness evolution similar to 
the maps in Figure 2. The movie is based on 142 outputs 
finely spaced between z = 0.66 and 0.26. Examination of 
the evolution showed that the large-scale cold front arises 
when a merging dense subclump is moving through hot- 
ter ambient intergalactic gas of the merging system. The 
hotter temperatures are due to the merger shocks^. The 
shocks are driven by supersonic motions of the merging 
clusters and thus propagate away from them (the situa- 
tion is similar to that of a supersonically moving piston 
in a gas tube). The shocks appear to originate as cluster 
pass each other for the first time when their velocity is 

-•^Bialek et al. (2002) pointed out that the gas sloshed out of its 
potential well can expand and adiabatically cool, which would also 
enhance the temperature difference between clump's and ambient 
gas. 
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Fig. 5. — Top panels: the map of projected dark matter density 
with overlaid contours of the emission-weighted gas temperature 
near the cold fronts identified in the z = 0.43 (left) and z = ^ 
(right) panels in Figure 2. The dark matter density maps are color- 
coded on a log^^Q scale in units of g cm~^; the temperature contours 
are spaced on linear scale. Bottom panels: the velocity maps of 
gas are overlayed on the emission-weighted temperature map of the 
regions shown in the top panels. The length of the thick vertical 
vector in the bottom-left corner corresponds to 1000 km s~^. The 
size of the shown region is 0.78/i~-'^Mpc. The velocity maps show 
that the cold fronts are non-equilibrium phenomena arising during 
cluster mergers. 

greatest. The clumps thus trail the bow shocks they pro- 
duce and the cold gas of clusters is not actually shocked 
until much later when the cold front is disrupted. Two 
additional weaker bow shocks are produced when the clus- 
ters pass each other for the second time. These secondary 
shocks play a role in heating the cold gas and disrupting 
the front. As the cold front observed in the temperature 
map is located at the large distance from the cluster cen- 
ter the front exists approximately for a time it takes for 
its sub-cluster to travel back to the center and merge plus 
the time required for the secondary shock to propagate 
from the center to the front. This is roughly equal to 
the dynamical time of the cluster (~ 10^ years). Kelvin- 
Helmholtz instability and merger-induced turbulence also 
appear to contribute to eventual disruption of the front. 
Based on this estimate, the frequency of large cold fronts 
should be somewhat smaller than the frequency of major 
mergers in progress because the cold front may appear only 
in certain projections (projection direction perpendicular 
to the plane of the merger). Cold fronts thus appear to 
be a common dynamic phenomenon accompanying merger 
events during cluster formation. 

4.2.2. Small-scale Cold front from Minor merger 

Examining outputs of the simulations, we find that 
large-scale sharp gradients of gas temperature most of- 
ten correspond to the boundaries between the relatively 
cold intracluster gas of merging clumps and hot, already 
shocked gas surrounding them. This is in agreement with 
interpretation of Markevitch et al. (2000) and Vikhlinin 



et al. (2001). We found many smaller cold fronts in fre- 
quent minor mergers. For example, a small-scale cold front 
associated with the merging subclump to the left of the 
cluster center can be seen in the z = panel of Figure 2. 
This cold front has a size of ^ bOh~^ kpc and temperature 
of ^ 2 keV, so this is essentially a poor group. Here, un- 
like the advanced merger stage discussed in the previous 
section, the merging subclump is on the first approach to 
the cluster center and has not yet suffered major tidal dis- 
ruption. The motion of the dense clump of gas compresses 
the ICM upstream of the front. The compression sharpens 
and enhances the amplitude of gas density and tempera- 
ture gradients across the front. The surface brightness, 
temperature and pressure profiles across this cold front 
are shown in the right panels of Figure 3. Here again, the 
simulated profiles reproduce all the main features of the 
observed cold fronts. 

The smaller fronts would be much more difficult to ob- 
serve because they have both much lower X-ray luminos- 
ity and surface brightness contrast; however, they can be 
observed in the nearby clusters. Recently, Forman et al. 
(2002) presented examples of small-scale cold fronts in the 
Fornax cluster which have a size of ~ 20 — 40 kpc and tem- 
peratures of < 1 keV. These observed properties are very 
similar to those found in our simulations. Note also that 
the relatively cooler gas of the merging sub-clump is trailed 
by low-entropy (~l-2 keV) intergalactic gas in the direc- 
tion of the subclump's motion (see also the entropy map in 
Figure 1). Such elongated regions of low-temperature gas 
trailing a subclump, if observed, can be identified with 
regions where filamentary material is accreting onto the 
cluster. The filamentary structures visible in the temper- 
ature map, however, does not appear in the X-ray surface 
brightness map. 

4.3. Dynamics of gas motions 

Detailed understanding of gas dynamics in the vicin- 
ity of cold fronts is important for placing meaningful con- 
straints on the physics of the ICM, since the power of these 
constraints relies on the validity of the underlying assump- 
tions of the dynamical model. 

Figure 5 presents projected DM density maps with over- 
layed contours of emission- weighted temperature (top pan- 
els) and emission-weighted temperature maps with over- 
layed gas velocity fields (bottom panels) for the two cold 
fronts discussed in the previous section. The velocity fields 
were constructed by averaging velocities in ^ 78h~^ kpc 
slice centered on the DM subclump associated with cold 
front. The figure reveals a variety of high- velocity streams 
arising during cluster merger. In the large cold front 
{z = 0.43) shown in the left panels, the dominant fea- 
ture of the velocity field is the collision of two gas fiows on 
the shock boundary. Note that there is no indication of a 
regular gas fiow around the cold front boundary. Instead, 
the cold gas of the front appears to be moving together 
with the surrounding hot gas. Moreover, it is clear that 
the subclump is in the process of tidal disruption: the ve- 
locity of its gas is highly non-uniform and the peak of the 
DM density is displaced with respect to the center of the 
front. In the z = minor merger (right panels) the ICM 
of the main cluster is considerably more relaxed (typical 
velocities of just ~ 100 — 300 km s~^). The bulk fiow 
of gas associated with the merging subclump is also quite 
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regular. Although in this case the assumption of regu- 
lar flow would be closer to truth, the figure shows that 
the flow on small-scales around the cold-hot interface is 
still quite complicated which may hamper development of 
Kelvin-Helmholtz instabilities. 

Finally, we would like to note that we performed a sim- 
ilar search for cold fronts in the SCDM cluster simulation. 
We found several cold fronts with properties similar to 
those found in the ACDM run. Here again the features 
were associated with mergers in different stages. In the 
two clusters that we analyzed, we identifled many "small- 
scale" cold fronts many of which are located in the out- 
skirts of the main cluster progenitor. 

5. DISCUSSION AND CONCLUSIONS 

We used high-resolution A^-body+gasdynamics cluster 
simulations in CDM models to search for counterparts of 
^'cold fronts within cluster ICM recently discovered in 
Chandra observations. The observed cold fronts are very 
sharp features in the X-ray surface brightness and temper- 
ature maps found in several clusters. These features are 
characterized by an increase in surface brightness by a fac- 
tor ^ 2 accompanied by a drop in temperature of a similar 
magnitude over the scale of 10-50 kpc. The existence of 
such sharp features can put interesting constraints on the 
properties of the intracluster medium (ICM) , including the 
efficiency of energy transport (Ettori & Fabian, 2000) and 
the magnetic field strength in the ICM (Vikhlinin et al., 
2001). In this paper, we investigate whether such features 
can be produced naturally in clusters forming in hierar- 
chical models and study their origin, structural properties 
and detectability. 

We find that cold fronts appear to be common and often 
arise in major and minor cluster mergers. In the preced- 
ing sections we discussed two specific cases of cold fronts 
found in the simulation of the ACDM cluster. In the most 
spectacular case, the cold front arises in a major merger 
when the merging subclump, moving slightly supersoni- 
cally, is undergoing tidal disruption (see distribution of 
DM and velocity field in Fig. 5). The motion of the gas, 
which appears to be in the process of being "sloshed out" 
of the potential well of its DM halo, drives the shock and 
leads to a steepening of density and temperature gradi- 
ents which produces the cold front. The surface bright- 
ness and temperature profiles across the front show that 
a drop of temperature over the scale of 20 — 30h~^ kpc is 
accompanied by a jump in surface brightness by a factor 
of ^ 2, behavior very similar to the profiles of observed 
cold fronts. In addition, we find that similarly to obser- 
vations pressure changes smoothly accross cold fronts in 
simulations. The spatial extent of the cold front in sim- 
ulation (^ 0.bh~^ Mpc) is comparable to the extent of 
two prominent observed cold fronts in clusters A2142 and 
A3667 (Markevitch et al., 2000; Vikhlinin et al., 2001). 
Note also that the cold front is preceded by a bow shock 
(~ 200 kpc from the cold front). In this case the jump 
in surface brightness by a factor of ^ 5 is accompanied by 
an increase in temperature of a similar magnitude over the 
scale of ~ 40 kpc. The shock front has a much lower X-ray 
surface brightness than the gas around the cold front and 
therefore would be difficult to observe. Nevertheless, in- 
dications of a similar bow shock preceding the cold front 
were found in A3667. 



In another discussed case, the cold front arises in a minor 
merger when a merging subclump reaches the inner regions 
of the larger cluster and gas in front of the subclump is 
significantly compressed. The compression sharpens and 
enhances the amplitude of gas density and temperature 
gradients across the front. Here the subcluster is on the 
first approach to the main cluster's core and is not yet 
significantly disturbed by tides. The relatively cooler gas 
of the merging sub-clump is trailed by low-entropy (^1- 
2 keV) intergalactic gas in the direction of the subclump's 
motion. Such elongated regions of low-temperature gas 
trailing a subclump, if observed, can therefore be identi- 
fied with regions where filamentary material is accreting 
onto cluster and could potentially be targets for studies 
of properties of intergalactic material. Unfortunately, the 
X-ray surface brightness of this gas is very low and de- 
tecting such regions will be a challenge. The best chance 
of observing such gas in candidate regions would be X-ray 
and UV metal absorption lines in spectra of bright back- 
ground sources (e.g., Hellsten et al., 1998). 

Cold fronts in both observed and simulated clusters 
seem to come in all shapes and sizes and are rather fre- 
quent, as reports of many new instances of cold fronts are 
published (e.g. Sun & Murray, 2002; Forman et al., 2002; 
Markevitch et al., 2002). In the two clusters that we an- 
alyzed, we identified many "small-scale" cold fronts many 
of which are located in the outskirts of the main cluster 
progenitor. Most of these cold fronts would probably not 
be detected in distant clusters with typical exposures and 
field sizes of Chandra observations, although they are sim- 
ilar to those observed in the Fornax cluster (Forman et al., 
2002). The fact that large-scale (- 0.5 Mpc) cold fronts 
arise in major mergers and survive for approximately dy- 
namical time (< 1 Gyr) implies that their frequency will 
be similar to the frequency of major mergers. The latter 
is a very strong function of redshift (oc (1 -h z)^^^ (see, 
e.g. Gottlober et al., 2001)) and may not be very high 
at z = (depending on cosmology and the power spec- 
trum normalization), but be increasingly larger for higher 
redshift clusters. Indeed, the large-scale cold front that 
we discussed occurs at z ^ 0.4. Observationally, cluster 
ZW3146 at z = 0.29 exhibits several cold front like fea- 
tures (Forman et al., 2002). More recently, Markevitch 
et al. (2002) reported results of a systematic search for 
cold fronts in archival Chandra observations of 30 appar- 
ently relaxed clusters. Surprisingly, they find cold- front 
like features in the large fraction of the clusters. Larger 
samples of both observed and simulated clusters will be 
useful to determine the frequency of cold fronts. 

Our results support the interpretation of cold fronts as 
boundaries between hot ICM of the main cluster and "the 
dense subcluster core that has survived a merger and ram 
pressure stripping by the surrounding shock-heated gas" 
(Markevitch et al., 2000). Based on this interpretation, 
Vikhlinin et al. (2001) proposed the detailed dynamical 
model of the cold front observed in cluster A3667 assuming 
a spherical dense gas subclump in a spherical dark matter 
potential moving through the uniform ICM. This model 
was used to place constraints on the magnetic field config- 
uration and strength and plasma conductivity in the vicin- 
ity of cold front. The conductivity can be constrained sim- 
ply by the sharpness of the observed front if a certain life- 
time is assumed for the front. The constraint on the mag- 
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netic field configuration, on the other hand, is obtained 
from the fact that no indication of Kelvin-Helmholtz in- 
stabilities is present in observations. Such instabilities are 
expected to develop quickly for a regular high- velocity flow 
of hot surrounding gas around the dense subclump along 
the cold front boundary. The power of these constraints 
relies on the validity of the underlying assumptions. Of 
these, the assumption of regular flow is particularly worri- 
some, given the complicated flow structure apparent from 
entropy and temperature maps shown in Figs. 1 and 2. 

Indeed, Figure 5 shows that this assumption fails in the 
case of cold fronts identified in our cluster simulations. 
Our simulations show that dynamics of gas around the 
simulated cold front is much more complicated than that 
assumed in Vikhlinin et al.'s model. In particular, the 
flow of gas around the front is not laminar and in general 
is not parallel to the front. Vikhlinin et al., on the other 
hand, assumed a laminar flow similar to flow of gas around 
a blunt body. In this case the gas flow lines are parallel 
to the front - the conditions ideal for development of the 
Kelvin-Helmholtz instabilities. The velocity field shown in 
Fig. 5 is very different and the conditions are probably not 
as conducive to the Kelvin-Helmholtz instability as those 
assumed in Vikhlinin et al.'s model. On the other hand, 
the cold front in A3667 that they model has a much more 
regular appearance than the simulated cold front shown 
in Fig. 4. But, such detailed comparison is difficult be- 
cause processes discussed by Vikhlinin et al. would be 
poorly resolved in our simulations, which cannot resolve 
well perturbations smaller than 20 kpc. It is also difficult 
to draw conclusions about how generic the differences are 
(i.e., whether the flow pattern assumed in Vikhlinin et al. 
can arise in simulated clusters). These questions will be 
addressed in future studies using considerably higher reso- 
lution simulations and a large sample of simulated clusters, 
respectively. 

To conclude, our results indicate that cold fronts are 
non-equilibrium transient phenomena arising in cluster 
mergers. They occur when the relatively cold gas of the 
merging subclump is moving fast with respect to the hot- 
ter gas of the main cluster. Both simulated cold fronts 
discussed in the paper would be detectable in a Chandra 
observation of ~ 50 ksec exposure, if the clusters were 
located at a moderately low redshift (e.g., z <^ 0.1, see 
Fig. 4). Cluster mergers (especially minor) are common 
and the fact that we were able to identify several cold 
fronts in the simulations of only two clusters means that 
we can expect many more cold fronts to be discovered in 
the near future (e.g.. Sun & Murray, 2002). If the origin 
and properties of these features are understood, they can 
provide powerful constraints on the frequency of cluster 
mergers and properties of the intracluster medium. The 
analysis presented here is a first step in studying the origin 
and properties of cold fronts. The presented results should 
be useful in assessing the validity of the assumptions in the 
detailed modeling of cold fronts and developing more so- 
phisticated models. 
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